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VO ' Abstract 

. We provide new existence and uniqueness results for the discrete- 

time Hamilton (DTH) equations of a symplectic-energy-momentum 
1 (SEM) integrator. In particular, we identify points in extended-phase 

. space where the DTH equations of SEM integration have no solution 

c| [ for arbitrarily small time steps. We use the nonlinear pendulum to 

Q-T illustrate the main ideas. 
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1 Background 

Is symplectic-energy-momentum integration well-posed? Loosely speaking, 
the answer is no. Points exist in the extended phase-space of a Hamiltonian 
system where the equations of a symplectic-energy-momentum (SEM) inte- 
grator have no solution for arbitrarily small time steps. Before considering 
this question in more detail, we provide a brief review of SEM integration. 

Hamiltonian dynamics is at the heart of modern physics and arises nat- 
urally in applications such as optimal control theory and geometric optics. 
Hamiltonian dynamics is also the inspiration for the relatively new field of 
symplectic geometry. A symplectic-energy-momentum (SEM) integrator is 
a numerical integrator that preserves the following key properties associated 
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with Hamiltonian dynamics: i) The integrator is symplectic. ii) The integra- 
tor exactly conserves energy (the Hamiltonian function), hi) The integrator 
exactly preserves "linear" symmetries (e.g. linear and angular momentum 
in Cartesian coordinates). The term "symplectic-energy-momentum integra- 
tor" was coined and popularized by Kane, Marsden and Ortiz [13]. See also 
Chen, Guo and Wu pQ for related work on higher-order, symplectic-energy 
integrators. Guibout and Bloch ^T] have developed a general framework 
for deriving many of the published symplectic integrators, including SEM 
integrators. 

The author's work on SEM integration — known as discrete-time, Hamil- 
tonian (DTH) dynamics — predates the work of Kane, et al. • DTH 
dynamics originated from an effort to obtain the exact energy and momen- 
tum conserving properties of the discrete mechanics of Greenspan [S], [5], 
from the variational principle used in the discrete mechanics of Lee [T5] > [E] • 
DTH dynamics was proved in 1994 (see Shibberu [12], [H]) to be symplectic 
and hence a SEM integrator. 

In the extended-phase space formulation of Hamiltonian dynamics, time 
is treated as a generalized coordinate on equal footing with the position 
coordinates. The momentum conjugate to time is introduced as an additional 
generalized coordinate. The principle of least (stationary) action takes a 
particularly simple form in extended-phase space. But, despite its aesthetic 
appeal, the extended-phase space formulation of the principle of least action 
is not widely used because it leads to indeterminate equations of motion [H] , 
[Hj, [2I]- Lee [T3], [15], described a discretization of Lagrangian dynamics that 
appeared to remove this indeterminacy. D'Innocenzo, Renna and Rotelli [3 
modified Lee's discretization and achieved exact energy conservation. SEM 
integration is based on a related, but more general, discretization developed 
independently of D'Innocenzo et al. [3] in Shibberu [TH] . 

An important theorem due to Ge (see [I] and citation in [E]) illustrates 
the difficulty of formulating a symplectic integrator which exactly conserves 
energy. Roughly speaking, Ge's Theorem says that a general, energy con- 
serving, symplectic discretization of Hamiltonian dynamics, must reproduce 
a reparametrization of the exact dynamics. Why SEM integration does not 
violate Ge's Theorem was explained for the first time in Shibberu [20] - 1 

This article is concerned with the following question. Under what condi- 
tions are the DTH equations of SEM integration well-posed? We will prove 
results which generalize the existence and uniqueness results first proved in 
Shibberu |18j . The existence and uniqueness results in this article are for 
nonlinear Hamiltonian systems and are local in nature. A global result for 
linear Hamiltonian systems was proved in Shibberu |18j . [2*T] . 

lr The explanation given in Kane et al. \1'6\ of why symplectic-energy-momentum in- 
tegration does not violate Ge's Theorem is incorrect. A variable-time step symplectic 
integrator can be reformulated in extended-phase space as a constant-time step symplec- 
tic integrator. Therefore, Ge's Theorem holds true even for variable time-step symplectic 
integrators. See the discussion of Hairer's J2] "meta-algorithm" for variable time-step 
symplectic integrators given in the last section of Shibberu [2] • 
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2 Example: The Nonlinear Pendulum 



In this section, we illustrate the main ideas of this article using the nonlinear 
pendulum as an example. We begin by describing how the DTH equations 
of Hamiltonian dynamics are derived. Then we consider the existence and 
uniqueness of solutions to the DTH equations. 

Let z = (q,p) T where q = (qi, . . . , q n , t) T and p = (pi, . . . ,p n , p) T are the 
extended phase space, position and momentum coordinates of an n degree- 
of-freedom Hamiltonian dynamical system with Hamiltonian function TC(z). 
The position coordinate t represents time and the momentum coordinate 
p represents the momentum conjugate to time. (See [T3j, or [T5] for 
a detailed description of p.) We represent the motion of a discrete-time 
Hamiltonian dynamical system by a piecewise-linear, continuous trajectory 
in extended-phase space where k = 0, . . . , N are the vertices of the tra- 
jectory and ~z~k, k — 0, . . . , N — 1 are the midpoints of the linear segments of 
the trajectory. 

Define the one-step action of a discrete-time Hamiltonian dynamical sys- 
tem to be the function A(zk, Zk+i) = §Ag fc T Ap fc . (The motivation for choos- 
ing this definition for the discrete action is given in Shibberu [22]-) The 
dynamics of a discrete-time Hamiltonian dynamical system is determined by 
the following variational principle. 

Definition 1 (DTH Principle of Stationary Action) The one-step ac- 
tion A(zk, -Zfc+i), k = 0, 1, . . .N — 1, is stationary along a DTH trajectory 
for variations which fix q^ and Pk+i and satisfy the Hamiltonian constraint 
H{z k ) = 0. 

The DTH equations of SEM integration are determined by Definition ^ 

Theorem 2 (DTH Equations) A DTH trajectory is determined by the 
following equations: 

Az k = \ k JH z {z k ) (la) 

H(z k ) = (lb) 



where J = 




and I is the n + 1 dimensional identity matrix. 



Theorem[2]is proved in Shibberu j22]- See also Shibberu [21] for the proof 
that the DTH equations ([la| ) -([lb |) preserve symplectic-energy-momentum 
properties and are coordinate invariant under linear symplectic coordinate 
transformations. 

For sufficiently small time steps, a sufficient condition for the existence 
and (local) uniqueness of solutions to equations ([la| ) -([TE ]) is the condition 
ip(zk) 7^ where ip = {J'^z) T 'Hzz{J'Hz) Shibberu [IB] . The new existence 
and uniqueness results proved in this article include points where ip(zk) may 
equal zero, but the Poisson bracket [ip, TC] \ z is not equal to zero. Smoothness 
requirements on the Hamiltonian function are also weakened from H. eC 3 (L r ) 
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Figure 1: A DTH trajectory of a nonlinear pendulum. The v-shaped curves 
correspond to points where ip(z) = an d the horizontal and vertical lines 
correspond to points where [ip,1~i\ = 0. 



to Ti. eC 2 (£7) where U C 9ft 2n+2 is on open set in extended-phase space. (See 
Theorem ITU on page ITU for the main result of this article). 

Consider now a nonlinear pendulum with extended-phase space Hamilto- 
nian function Tt(q,p, p) = p + \p — cos(g). (Recall that p is the momentum 
conjugate to time.) The corresponding discrete-time Hamilton (DTH) equa- 
tions are 



Figure UJ is a plot of a DTH trajectory determined by the above equations 
and projected onto the phase portrait of the pendulum. Observe that the 
linear segments of the DTH trajectory are tangent to an energy conserving 
manifold of the pendulum. (We stress that the size of the initial time step, 
Ao, is determined by the initial condition zq = (qo,to,Po,po)-) The v-shaped 
curves in Figure ^ are points where if)(z) equals zero. The horizontal and 
vertical lines are points where the Poisson bracket [if>,H] equals zero. From 
Figure UJ we see that the existence and uniqueness results in this article 
apply to all the points in phase space except the equilibrium points where 
both ip(z) and [ip.'H] are equal to zero. 

Let ipk = tp(zk)- We will show that, for points where if>k 7^ 0, the mag- 
nitude and sign of Tik/^k is key to determining if a solution to the DTH 
equations exists and is locally unique. In particular, if TLk/ipk < 0, and 
ipk is sufficiently large, then no solution exists. If ipk = 0, the quantity 
Tik/ipki where ip' = [ip,7{\, plays a similar role in determining existence and 
uniqueness. In the neighborhood of points where if) changes sign, a DTH 



At k 



A*, 

-A fe sin(g fc ) 




Pk + 7^1 ~ cos(q k ) 



0. 



trajectory bifurcates giving rise to "ghost trajectories". Ghost trajectories 
are discussed in more detail in section [Tj 

The outline of this article is as follows. In section ??, we use the Newton- 
Kantorovich Theorem to prove the existence and uniqueness of a function 
z(X,z k ) implicitly defined by equation ((Taj). We use the function z(X,z k ) to 
decouple equation (JTb"]l from equation (JTa|). In section HI we derive a cubic 
approximation of the Hamiltonian constraint function g(X, z k ) = H(z(X, z k )). 
In sectional we identify intervals where g(X, z k ) is monotonic increasing/ decreasing 
with respect to A. Using monotonicity and the Intermediate Value Theorem, 
we prove the existence and uniqueness of Lagrange multipliers satisfying the 
decoupled, Hamiltonian constraint equation g(X,z k ) = 0. The existence and 
uniqueness results for Lagrange multipliers is used in section 13 to prove the 
existence and uniqueness of DTH trajectories. SEM integration is shown, 
under certain conditions, to be well-posed. Finally, in section (7j we discuss 
ghost trajectories and the need to regularize the DTH equations of SEM 
integration. 

3 Existence of a Decoupling Function 

Consider the DTH equations (jla|) - (jlbj) . Equation (jla|) can be rewritten 
as f(X,z k ,z k ) = z k - z k - \XJH z (z) = where z k = \{z k+1 + z k ). In 
Theorem below, we prove that if the Hamiltonian function 7i{z) satisfies 
certain conditions, then there exists a smooth function z(X, z k ) such that 
/(A, z k , z(X, z k )) = for all A G [—Xs, Xs] and z k G Us where Xs and Us are 
specified in Theorem El The function z(X, z k ) is used in section EJto decouple 
equation (jTb|) from equation ((Taj). We begin by stating two standard results 
in numerical analysis, the Newton-Kantorovich Theorem and the Matrix 
Perturbation Lemma [7j. 

Theorem 3 (Newton-Kantorovich Theorem) Consider the function f : 
U — > lZ n where U C 1Z n is open. Assume f G and \\f x (x 2 ) — fx(%i)\\ < 

7 \\x2 — Xi\\ for all X\,X2 G U. Assume there exists a point xq G U and con- 
stants (3 > 0, r] > such that \\fxixo) 1 \\ < P ^nd ||/ 3; (^o)~ 1 /( ;:c o) II < V- 
Assume a < | where a = (3^r]. Define r_ = (l — y/\ — 2a) / (3^ and r + = 
(l + y/1 — 2a) / Pj. If the close ball B(xo,r-) C U, then the Newton iterates 
x^\ defined by x^ 1 ^ = — f x {x^)~ 1 f{x^), i — 0, 1, . . ., with x^ = xo, are 
well defined and converge to x* G B(x ,r_) where is the unique solution 
of f{x) = in B(x , r + ) n U 

Lemma 4 (Matrix Perturbation Lemma) Assume the identity matrix 
I is perturbed by the matrix E. If\\E\\ < 1, then (I — E)~ exists, (I — E)~ = 
J2Z E" and ||(7 - E)' 1 }] < 1/(1-||S||). 

Theorem 5 (Decoupling Function) Consider the extended-phase space Hamil- 
tonian function H G C 2 {U) where U C TZ 2n+2 is open. Assume \\H Z {. Z )\\ < 
Mi and \\H zz {z)\\ < M 2 for all z G U. Assume \H zz {z x ) - H zz {z 2 )\ < 
1h \\z x - z 2 \\ for allz x z 2 G U Let X s = min(l/M 2 , 1/ 7h , (1 - (1 - 5) 2 ) j2M x ) 
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{z : B(z, S) C U} . Define /(A, z,z) = z — z — \XJTL Z (z) . Then, 



and Us 

there exists a 5, where < 5 < 1 
tiable function z : [—Xs, Xs] x Us - 
a// (A, 2) G [-Xg, X s ] x Efj. 



and there exist a continuously differen- 
TZ 2n+2 , such that f(X,z,z(X,z)) = for 



Proof. First we show that for |A| < 1/M 2 , f z 1 exists and is bounded. Since 
f z = I — E, where E = \X3H ZZ , and since \\E\\ < ± (1/M 2 ) M 2 = | < 1, 
by the Matrix Perturbation Lemma, J^ 1 exists and If/^ 1 ]) < P, where j3 = 
1/ (l — |) = 2. Next, we show that for |A| < 1/jh, h is Lipschitz with 
respect to z with Lipschitz constant 7 



~. We have 



||/»(A,z 2 ,«) - /»(A,zi,z)|| 



11 „_ _ 

< o lH \\Z2 ~ Zi\ 



\z 2 - Z\\ 



Now consider using Newton's iteration to solve /(A, Zk, z) = for z. If we ini- 
tialize the iteration with ~z^ = Zk, we have 77 
Wf-'W UXJH Z \\ . Let A 5 = min(l/M 2 , 1/ 7h , 



fz (X, Z k , Zk)f(X, Zk, z k ) < 



(1 



/2Afi). For |A| < 

X s , it follows that 77 < 2|(1 — (1 — 5) 2 ) / (2M 1 ) M x = (1 - (1 - 5) 2 ) /2. For 
< 5 < 1, we have then that a = /3 7 ?? < 2| (1 — (1 — 5) 2 ) /2 < ± Therefore, 



<5) 2 ) 
fl- 



1 - VI -2a 
r_ = < 

Pi 



1 - Wl 



i-(i-g) 2 



2± 



It follows that for z k G t/j, B{zk,rJ) C B(zk,S) C [/. By the Newton- 
Kantorovich Theorem, the function z(X, z k ) is well defined on [—Xs, Xs] x C/^ 
and /(A, ^(A, ^fc)) = 0. (We assume 5 is chosen small enough that U$ is 
nonempty.) The Implicit Function Theorem implies z(X,Zk) is continuously 
different iable. ■ 



4 Cubic Approximation of the Hamiltonian 
Constraint 

Given X k G [—Xs, Xs] and G Us, Theorem El implies there exists a point 
Zk = z(Xk, z^ and a point Zk+i = %Zk— Zk, such that Afe, Zk and Zk+i satisfy 
the first DTH equation, Azk = XkJTC z (z~k) . We use z(X, Zk) to decouple the 
second DTH equation, TC(zk) = 0j from the first DTH equation by defining 
the function g(X,Zk) = TC(z(X, Zk)) and replacing the second equation with 
the equation g(A, Zk) = 0. 

In this section, we determine a cubic approximation of g(X, Zk) as a func- 
tion of A. Obtaining this approximation is made difficult by the fact that the 
function z(X, Zk) is only implicitly defined. We will see that the linear term 
in the cubic approximation of g(X, Zk) is always equal to zero. The analysis of 
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DTH dynamics is also complicated by this fact since we are forced to consider 
the effects of the quadratic and even cubic term in the cubic approximation 
of g{X,z k ). 

The outline of this section is as follows. In Lemma |B1 below, we show that 
z\(X,z k ) is Lipschitz continuous with respect to A. In Lemma d we define 
the important function ip(z) = (JH Z ) T H zz (JH Z ) and we approximate the 
partial derivative dg(X, Zk)/d\ by the simpler function — jXh(X, z k ) where 
h(X, Zk) = ip(z(X, z k )). In Lemma |Hl we prove that dh(X,z k )/dX is Lipschitz 
continuous with respect to A. Finally, in Lemma El we determine a cubic 
approximate of g(X, z k ). 

Lemma 6 For Ai, A2 G [— Xs, \g] and z k G Us, 

\\z\(\2,Zk) -z x (Xi,z k )\\ < j z |A 2 - Ai| 
where j z = 2M X M 2 + M 2 . 
The proof is given in the appendix. 

Lemma 7 Define the functions g(X, z k ) = H(z(X, z k )), if)(z) = (JH Z ) T H zz {JH. Z ) and 
h(X,z k ) = i[)(z(X, z k )). Then, for |A| < Xs and z k G Us, 

< -M?Ml |A| 3 . 
- g 1 211 

Proof. Since /f 1 = (/ — E)~ l where E = \XJH ZZ , for |A| < X 5 we 
have \\E\\ < | < 1. By the Matrix Perturbation Lemma we have 

j= 1 = / + E + E 2 + E 3 + E A + ■ ■ ■ 
= I + E + E 2 + E 3 (I + E + ■ ■ ■ ) 
= I + E + E 2 + E 3 f-\ 

Therefore, 

dg(X, z k ) _ T _ 

~^x~- HzZx 

= \nj (f-^J) n z 

= \{H]JH Z + H] (EJ)H Z 

+ H] (E 2 J) H z + H] (E 3 fs\j) H z ) . (3) 

Since both J and E 2 J = ^X 2 (JTi. zz JTi, zz J) are skew-symmetric, the first and 
third term in (j3J) equal zero. The second term is given by 

H T Z (EJ) n z = - l -X (JH Z ) T H zz (JH Z ) 

= -o^V^OMfc)) 

= - l -Xh{X,z k ). (4) 



dg{X,z k ) 
dX 



-~Xh(X,z k ] 
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Thus, equations © and 
dg(X, z k ) 



OX 



imply 



--A/i(A, z k ) 



l -Hl(E*f^j) H z 



< -M\Ml |A| 3 . 



Lemma 8 Assumeip £ C 2 (U), \\ip z (z)\\ < Ni and 
Then, for Ai, A 2 £ [—\s, Xs] an d Zk £ Us, 



[z)\\ < N 2 forze U. 



dh(\ 2 ,z k ) dh{\ u z k ) 



< lh |A 2 - Ai 



d\ OX 
where j h = N^ z + M?N 2 . 
The proof is given in the appendix. 



Lemma 9 Assume ip £ C 2 {U), \\ijj z {z)\\ < N± and \\ij) zz {z)\\ < ^"2 for z £ U 
Let Hk = H(zk), ipk = ip(zk) and ip' k = [ip,H)\ z=z ■ Then, for |A| < Xs and 
Zk £ U s , 



dg(X,z k ) 



and 



dX 



g{x,z k ) 



-l^x-^kX 2 



n k - -i>kX 2 - ^' k x 3 



1 

24 



< AK\ 



<K\ 



(5) 
(6) 



where K = ^(MfMf + 2 7h ). 



Proof. The Mean Value Theorem implies there exists a A between and A 
such that h(X, z k ) — /i(0, z k ) = (dh(X, Zk)/dX) X. Therefore, using Lemma |H1 



dh 

h(X, z k ) - h(0, z k ) - ^(0) z k) A 



dh ~ . dh 

— {X,z k )X- —{0,z k )X 



<i h \x\\x\< lh \xf 



Since h(0, Zk) = ipk and dh(0, Zk)/dX = \ip' kl we have \h(X, z k ) — ipk — \^' k A < 
\ X\ 2 ■ Using Lemma [7[ we establish inequality (JSJ) as follows. 



dg(X,z k ) 



OX 



~^kX- -i>' k X 2 



< 



+ 



dg(X,z k ) 



OX 



-~Xh(X,z k ) 



-~Xh(X,z k ) 



-\^kX-^ip' k X 2 



< ^K\X\ : 



(7) 



S 



where K = -UMfMf + 2ry h ). We establish © as follows. First, using © 



we have 

A 



dg(X,z k ) 1 



I \2 



dX 4 



+ + -v4 a 



rfA 



<AK \Xr dX 



K\X\ 



But 

»A 



J (^%^ + \^ X + l^ X2 ) dX = g(X,z k )-g(0,z k ) + ^ k X 2 + ^' k 



X". 



So we have 



<i(X. :/,) - ( 7Y fc - i^A 2 - l^A 3 



5 Existence and Uniqueness of Lagrange Mul- 
tipliers 

In this section, we address the question of the existence and uniqueness of 
Lagrange multipliers A which satisfy the decoupled, Hamiltonian constraint 
equation g(X, z k ) = 0. We begin by proving a monotonicity result for the 
function g(X,z k ). Then we prove three separate existence and uniqueness 
theorems, Theorems ITTHT^l each of which accounts for one of the three regions 
of extended-phase space described below. (The value of the constant K is 
determined by the Hamiltonian function H(z). See Lemma 03) 

region I {z k : ijj{z k ) ^ 0, {ip'{z k )f < 2AK \ip{z k )\] 
region II {z k : ${z k ) + 0, (f(z k )) 2 > 24K \ip(z k )\} 
region III {z k : ip(z k ) = 0, ip'{z k ) ^ 0} . 

The proofs of each of the three existence and uniqueness theorems in 
this section uses the same basic approach. First, we derive bounds for the 
function g(X, z k )/ip(z k ). (See Figured) Then, we use monotonicity and the 
Intermediate Value Theorem to establish the existence and (local) uniqueness 
of Lagrange multipliers A satisfying the equation g(X, z k ) = 0. 

Lemma 10 (Monotonicity) Assume z k G Us- Then we claim the fol- 
lowing: 

(%) If fa 0, (^) 2 < 2AK \fa\ and |A| < A k where < A k < min(^/\fa\ /96K, X s ), 
then g(X, z k ) is monotonic increasing/decreasing in the intervals (—A k , 0) 
and (0, Afe). 
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Figure 2: Bounds on g(X, z k ) /ip(z k ) for the nonlinear pendulum. 




Figure 3: Plots of g(X, z k )/ip(z k ) in region I of the nonlinear pendulum. 



(ii) Assume ip k ^ 0, (ip' k ) 2 > 24K \ip k \ and \X\ < A k where < A k < 
mm(\ip' k \ /48K, Xs). Let g(s, z k ) be a reparametrization of g(X, z k ) where 
s = — (ip' k /ip k )X. Define S k = \ip' k /ip k \A k . Then g(s,z k ) is monotonic 
increasing/decreasing in the following intervals: a) (—S k ,0), b) (0,S k ) 
ifS k < I c) (0,f) ifS k > | and d) (6, S k ) if S k > 6. 

(Hi) Ifip k = 0, ip' k ^ 0, and \X\ < A k where < A k < mm(\ijj' k \ /4SK, Xs), 
then g(X, z k ) is monotonic increasing /decreasing in the intervals (—A k , 0) 
and (0, Afc). 

The proof of Lemma El is given in the appendix. 

Theorem ^2 below deals with region I of extended-phase space where 
the quadratic term dominates the cubic term in the cubic approximation of 
g(X, z k ). See Figure El for plots of g(X, z k )/ip k in region I of the nonlinear 
pendulum. Since g(0,z k ) = H(z(0, z k )) = TL(z k ) = Ti k , we see from FigureEl 
that the sign of TLk/ipk determines the number of solutions to the equation 
g{X,z k ) = 0. 

Theorem 11 Assume z k e Us, ip k ^ 0, (ip' k ) 2 < 24K \ip k \ and \X\ < A k 
where < A k < min(^\ilj k \ /9QK, Xs). Then the following statements about 
the equation g(X, z k ) = are true. 
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(i) IfH k /ip k < 0, no solution exists. 

(ii) IfHk/ipk — 0, the only solution is X = 0. 

(Hi) If < H k /ip k < JjjAjfc, two solutions of opposite sign exist, X k G 
(— Afc, 0) and X k G (0, Ajt). T/ie solutions are unique within their re- 
spective intervals. 

(iv) IfH k /ip k > ;i>A^, no solution exists. 



Proof. Since |A| < A k < 
EJthat 



/96K, we have from inequality (jBJ) of Lemma 



H k 



*pk 

It follows that 

V± _ I^a 2 - — ^ 
^ 96 24 ^ 

Since by assumption (ip' k 



_ X 2 _ 1^3 

S 24^ 



< 



|A| 4 <-A 2 
1 1 ~ 96 



A 



<24K\^ k 



X 2 < flCMfc) < ?4 

we have 



1. 2 1 

— A H 

96 24 



^k 



A 



^A 













< 



96K 



1 

2' 



Using (JHJ) and © we have 



V»fe 32 



< 



g{\z k ) 



< 7j± 3^^ 2 



(8) 



(9) 



(10) 



ijj k ifj k 32 

To establish (i), assume H k /ipk < 0. Then inequality (fTU|l implies g(X, z k )/ip k < 
for all |A| < A k and no solution exists. If H k /ip k = 0, then g(X, z k )/ip k < 
for nonzero A. Since g(0, z k ) = H k = 0, the only solution is A = 0, establishing 
(ii). If we assume < H k /ip k < ^A 2 ,, then g(±A k , z k )/tp k < H k /ip k - ^A 2 k < 
0. Since g(0,z k )/ip k = H k /ip k > 0, the Intermediate Value Theorem implies 
g(X,z k ) = has two solutions X k G (— A fc ,0) and X k G (0,A fc ). Lemma ITUlT j 
implies g(X, z k ) is monotonic in each interval establishing uniqueness and 
claim (iii). Finally, if H k /ip k > J^A 2 ,, then inequality (fTUj) implies that for all 
IAI < A*. 



< 



Hi 



-A . < 



Hi 



-A 2 < 



^ k 32 K ~ ^ k 32 " ^, k 
establishing claim (iv). ■ 

Theorem ^] below deals with region II of extended-phase space where 
i/j(z) is small but nonzero. Theorem is the most complex of the three 
existence and uniqueness theorems in this section because both quadratic 
and cubic terms need to be taken into consideration. The reparametrization 
s = — (ip' k /ipk) X simplifies the statement of the theorem and its proof. See 
Figure 13] for plots of g(X, z k )/ijj k in region II of the nonlinear pendulum. We 
can see from Figure 0] that the sign of H k /ip k determines the number of 
solutions to the equation g(s, z k ) = 0. (Recall that o(0, z k ) = H k .) Note the 
appearance of a "ghost solution", s* k . From Figure |4(cj| we see that, unlike 
the two solutions s k and s k , the ghost solution si does not approach zero as 
H k /v k -+ ()-. 
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Figure 4: Plots of g(X, Zk)/ip{zk) in region II of the nonlinear pendulum. 

Theorem 12 Assume z k G Us, ipk ^ 0, (ip' k ) 2 > 2AK \i/j k \ and \X\ < A k 
where < A k < min(\ip' k \ /48K, Xs). Let g{s,z k ) be a reparametrization of 
g(X, z k ) where s = — (ip k /ip k ) X. Define S k = \ip' k /ip k \ A k . Then the following 
statements about the equation g(s, z k ) = are true. 

(i) IfH k /ip k < 0, no solution exists in the interval (—S k ,2). 

(ii) IfH k /ip k = 0, the only solution in the interval (—S k ,2) is s = 0. 

(Hi) If ^A^(6 — S k ) < H k /ipk < for S k > 6, there exists a solution 
s* k e[2,S k ). 

(iv) IfO < TCk/ipk < 4^Afc(6 + Sk), there exists a solution s k G (—S k ,0) and 
the solution is unique in this interval. 

(v) If < TCk/ipk < 76^(2 — Sk) for S k < |, there exists a solution 
s k G (0, S k ) and the solution is unique in this interval. 

(vi) IfO < Hk/^k < ifg (^M) 2 and 

(a) if Sk > §, there exists a solution s k G [0, |) and the solution is 
unique in this interval. 

(b) if Sk > 6, there exists a solution s* k G 

(vii) IfH k /ip k > § (V'fc/Vi) 2 ? no solution exists in (0,S k ). 

The proof of Theorem ^] is given in the appendix. 

Theorem EH below deals with region III of extended-phase space where 
the quadratic term of the cubic approximation of g(X, z k ) is equal to zero. See 
Figure for plots of g(X, z k )/ip' k in region III of the nonlinear pendulum. As 
we can see from Figure the sign of 7i k /ip' k determines whether the solution 
A^ or A^ exists. 

Theorem 13 Assume Zk G Us, ipk — 0, ip' k 7^ an ^ \M — ^-fc where < 
Ak < mind^l /48K, Xs). Then the following statements about the equation 
g(X, z k ) = are true. 
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Figure 5: Plots of g(X, Zk)/ipf(zk) in region III of the nonlinear pendulum. 

(i) If — < Tik/ip'k < 0, tnere exists a solution A fc G (— Afc,0) and it is 
unique in this interval. No solution exists in [0, A&). 

(ii) If < TCk/ip k < ^A^, tnere exists a solution A^ G (0, A&) and it is 
unique in this interval. No solution exists in (— Ajt,0]. 

(Hi) If Tiki Vfc = 0; the only solution is A = 0. 

fiu) If\Hk/4>' k \ > i^A^, no solution exists. 

The proof of Theorem EI] is given in the appendix. 



6 Existence and Uniqueness of DTH Trajec- 
tories 

The main result of this article is stated below in Theorem The proof 
of Theorem [21 uses Theorems HTUT31 from the previous section. Before stat- 
ing the theorem, we provide a condensed description of the theorem's main 
conclusions. 

Consider a point z in extended-phase space. Roughly speaking, when 
|7io/Vo| is sufficiently small, there are four generic possibilities for DTH tra- 
jectories. (1) If Ho/^o > and \ipo\ is large, then a unique DTH trajectory 
exists which passes through the vertex point Zo. (2) If TCo/ipo > and if) 
changes sign near zq, then a DTH trajectory exists which bifurcates at the 
vertex point zq. (3) If Ho/?po < and ip changes sign near z Q , then a DTH 
trajectory exists which either begins or ends at the vertex point zq. (4) If 
H-o/ipo < and \ipo\ is large, then no DTH trajectory can exist having z as a 
vertex point. See FigureElfor plots of DTH trajectories of the nonlinear pen- 
dulum for the following initial conditions: (a) zq = (—0.15, —0.05,0, —3.48) 
(b) z = (1.25, -0.05, 0, -4.33) (c) z = (1.85, -0.05, 0, -4.87). 

Since DTH trajectories preserve the symplectic-energy-momentum prop- 
erties of Hamiltonian dynamics, Theorem IT41 provides conditions under which 
a SEM integrator is well-posed. As a practical matter, we point out that, 
for classical Hamiltonians, generic possibility (4), where no DTH trajectory 
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Figure 6: DTH trajectories of the nonlinear pendulum. 

exists, can always be avoided by choosing an initial value for po (the mo- 
mentum conjugate to time) which is sufficiently small and of the appropriate 
sign. Generic possibilities (2) and (3) are more challenging to deal with and 
are discussed further in section [7| 

Theorem 14 (Existence & Uniqueness of DTH Trajectories) Consider 
an extended-phase space Hamiltonian function 7i E C 2 (U) where U C 3? 2n+2 
is open. Define ip{z) = (JH Z ) T H zz (JH Z ) and ip'{z) = [ip,H}. Assume H z , 
T~Lzz, ifrz an d ip zz are bounded on U and TC ZZ is Lipschitz continuous on U. 
Assume also that Zk E Us where Us — {z : B(z, 5) C U\ and < 5 < 1. // 
ipk 7^ 0, then there exists a > for which statements (i)-(iii) are true. 

(i) IfTi,k/ipk < 0, then Zk can not be a vertex point or end point of a DTH 
trajectory with Lagrange multiplier (s) |A| < 

(ii) IfTik/ipk = 0, then Zk is a vertex point of a fixed-point DTH trajec- 
tory with Lagrange multiplier X k = 0. No other DTH trajectory with 
Lagrange multiplier (s) |A/-| < A k exists. 

(Hi) IfTik/ipk > is sufficiently small, then Zk is a vertex point of a unique 
DTH trajectory passing through Zk with Lagrange multipliers |A^| < A^. 

If IV'&I 7^ sufficiently small and ip' k ^ 0, then there exists a Ak > for 
which statements (iv)-(vi) are true. 

(iv) If Tik/ipk < and \7ik/ipk\ is sufficiently small, then Zk is a vertex 
point of a DTH trajectory which begins (ends) at Zk and has Lagrange 
multiplier \\* k \ < Ak- 

(v) IfHk/ipk = 0, then a DTH trajectory exists which bifurcates at Zk into 
a, fixed point, DTH trajectory with Lagrange multiplier Xk = and a 
ghost DTH trajectory with Lagrange multiplier |A£| < A^. 

(vi) IfHkj^k > is sufficiently small, then a DTH trajectory exists which 
bifurcates at Zk into a DTH trajectory with Lagrange multipliers | A^| < 
Ak and a ghost DTH trajectory with Lagrange multiplier |A£| < A&. 
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If 4>k = and ip'k 7^ 0, then there exists a A k > for which statements (vii) 
and (viii) are true. 



(vii) If \H k /ip' k \ 7^ is sufficiently small, then z k is a vertex point of a 
unique DTH trajectory which begins (ends) at z k with Lagrange multi- 
plier \X k \ < Aft. 

(viii) IfTC k /ip k = 0, then z k is a vertex point of a fixed-point DTH trajec- 
tory with Lagrange multiplier X k = 0. No other DTH trajectory with 
Lagrange multiplier |Aft| < A k exists. 

Proof. Consider the DTH equations 

Az k = XJH z {z) (11) 
H{z) = (12) 

where Az k = z — z k , z = \{z + z k ) and H EC 2 {U). We can rewrite equation 
(fTTj) as follows: 

f(k,z k ,z) = z- z k - -XJH z {z) = 0. (13) 

By assumption, there exists Mi, M 2 and 7# such that ||W Z (^)|| < Mi and 
\\H zz (z)\\ < M 2 for z E U and ||W„(zi) - H zz {z 2 )\\ < 1h\\z 1 -z 2 \\ for 
Zi,z 2 E U. Define 

Aj = mm — , — , 



•U, -„ 2Mi 7 

Theorem El implies there exists a < 5 < 1 and a function z(X, z k ) such that 
/(A, z k ,z(X, z k )) = for all (A, z k ) E [— X$, kg] x L^. Use the function z(k, z k ) 
to decouple equation (|T2jl from equation (jTTj) to obtain the equation 

g(X,z k )=H(z{X,z k )) = 0. (14) 

For a given z k E Us, equation (jTjj) determines the value of a Lagrange mul- 
tiplier Aft, provided that one exists. If a Lagrange multiplier(s) exists, then 
Aft and z k determine z k _i and/or z k+ i as follows: 

If Aft < 0, define z fc _i = z(X k , z k ) and z fc _i = 2z k -i - z k . 
If Aft > 0, define z k = z(X k , z k ) and z k+1 = 2z k - z k . 

The extended-phase space, vertex points z k -x, z k and z k+ i, determine a DTH 
trajectory which passes through the vertex point z k . If only z k _\ or z k+ i 
exists, then the DTH trajectory either begins or ends at z k . 

Now, we consider the existence and uniqueness of solutions to equation 
(|Hjl . By assumption, there exists N\ and N 2 such that ||^ z (2r)|| < N\ and 
Uzz{z)\\ < N 2 for zEU. Define 

K = — (Mf M| + 1M.U. Y ; + 2M 1 2 N 1 + 2MlN 2 ). 
32 
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Assume ip k ^ 0. If (ip' k ) 2 < 2AK \^ k \ , choose < A a < min(^|^ fc | /967T, X s ). 
Then Theorem ITT1 (i)-(iii) imply statements (i)-(iii) are true. If, on the other 
hand, (tfj' k ) 2 > UK\ip k \, choose < A b < min(\ip' k \ /48K, X 6: § l^fcMJ)- 
For this choice of A 6 , we have S k = \ip' k /ip k \ A b < |V4Mfc| § iV'fc/Vil = §■ 
Since S k < §, Theorem IT21 (i). (ii), (iv) and (v) imply statements (i)-(iii). 
Therefore, for A k = min(A a , A&), statements (i)-(iii) are true. 

Next, we prove (iv)-(vi). Assume < \ip k \ < min((^) 2 /288i^, \ip' k \ A5/6). 
Then 6 \ip k /ip k \ < min(|^| /A8K,Xs). Therefore, we can choose A k so that 
GhfoMI < A fc < min(|Vfe| /48K,X 5 ). Then we have S k = \ip' k /ip k \A k > 
Wkl^kl 6 iV'fc/V'fcl = 6. If 7i k /ip k < 0, Theorem W2X i) implies no solution exists 
in (— S k ,2). Since S k > 6, Theorem l^f iii) implies that for sufficiently small 
|Wfc/V'fc| , a solution s* k G [2,S k ) exists. Since A£ = —(ip k /ip' k )sl, and si > 0, 
all solutions that may exist must have the same sign. Hence, the DTH 
trajectory must begin (end) at z k , proving statement (iv). If H k /ip k = 0, 
Theorem lT2T ii). (vi)b prove statement (v). If Ti k /ip k > is sufficiently small, 
Theorem lT2T iv) . (vi) prove statement (vi). 

Finally, assume ip k = and ip' k 7^ 0. Choose < A k < mm(\ip' k \ /4SK, Xs). 
Theorem ^fi)-(iii), imply (vii)-(viii). ■ 

7 Ghost Trajectories 

Discrete approximations of differential equations can introduce spurious or 
nonphysical solutions. Greenspan (TH] provided a detailed analysis of a "non- 
physical" solution to his equations for discrete mechanics. Greenspan showed 
that, unlike the correct physical solution, the nonphysical solution approaches 
infinity as the time step is brought to zero. 

Multiple solutions also exist in DTH dynamics. When IV^)! is large, 
the decoupled Hamiltonian constraint equation g(X,z k ) = has only two 
solutions, X k and A^ , both of which appear to represent the correct physical 
behavior of the system. The Lagrange multiplier A^ corresponds to the 
trajectory propagating backward in time from z k and A^ corresponds to the 
trajectory propagating forward in time from z k . 

Near points where ip(z) changes sign, a third solution to g(X,z k ) = 
appears — the solution X k . As stated earlier, the solution A£ has a property 
that distinguishes it from the solutions X k and X k . Assume a sequence of 
initial conditions z k approaches the Hamiltonian conserving manifold T~i{z) = 
0, but not the manifold if)(z) = 0. Then the corresponding sequences, X k and 
X k , each converge to zero, but the sequence X* k does not converge to zero. 
We make this property precise in Theorem 1151 below. 

The solution X* k causes a DTH trajectory to bifurcate at z k , giving rise to 
what we call "ghost" DTH trajectories. Ghost DTH trajectories are not time 
reversible. (See Shibberu [22| for the details.) We will refrain from calling 
ghost trajectories "nonphysical" because, in DTH dynamics, it is unclear 
what the physically correct solution across ip(z) = manifolds should be. 
It appears that DTH dynamics needs to be regularized in some fashion. In 
Shibberu |22] > we propose a regularization of DTH dynamics which preserves 
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symplectic-energy- momentum properties and time reversibility across ip{z) = 
manifolds. 



Theorem 15 Consider a sequence Zk G Us, k = 0, 1,2, . . . where \ip(zk)\ > 
VVin > and TC(zk)/ip(zk) — > + . // A^ exzsi, i/ien lim^oo A^ = 0. // 
lim^oo X* k exists, then lim*^^ \\* k \ > § (VWn/MiiVx) > 0. 

Proof. Assume Zk is in region I of extended-phase space. Then, inequality 
of Theorem E] implies 

°4(^lSf (15) 

Now assume Zf. is in region II. Depending on the sign of ip' k /ipk, either s k = 
— (i)'k/if)k,) or s k = ~ O'k/i'k,) ^t- Therefore, inequality (|2~5j) of Theorem 
[T21 implies either 

Q <^(K) 2 (Q + h\)<^H (16) 

48 v k> v 1 fc u ip(z k ) 

or 

o^(*)'(e + |*|)<$g (17) 

Likewise, since s k < | < 2, inequality (|2T)j) implies, in correspondence with 
inequalities (fTH |) - (fT7]) . either 

0<i(A l -r ( 2- sJ )<^M (!8) 

or 

< ^ M)' (2 - s t ) < *g. (19) 

Since Hk/ipk — > + , inequalities (fT3j) - (fTT?j) imply lim^oo A^ = 0. Since |s£| = 
H^MKI > §, we have \\*\ > § |V*MI > f (VwM^i) . If lim^ A£ 
exists, we must have lim^oo |A£| > | (VWi/^iM) > 0. ■ 



8 Conclusions 

The extended-phase space formulation of the principle of least action leads 
to indeterminate equations of motion. Since SEM integration is based on 
a discrete version of this principle, it is important to establish conditions 
under which the equations of SEM integration are well-posed. Theorem ITU 
provides such conditions. Theorem ITU also shows that the DTH equations of 
SEM integration need to be regularized in some fashion. One proposal for 
regularizing SEM integration is given in Shibberu [2*2"] . 

The existence and uniqueness results in this article are only locally valid. 
A global result — for example, sufficient conditions for the existence of DTH 
trajectories for arbitrarily long intervals of time — would be interesting. One 
of the difficulties in establishing such a result appears to be establishing a 
global bound on the Lagrange multipliers A&, k = 0, 1, 
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A coordinate-invariant, formulation of DTH dynamics could provide ad- 
ditional insight into the behavior of DTH trajectories crossing ip(z) = 
manifolds. Preliminary work on a coordinate-invariant formulation of DTH 
was given in Shibberu ^H] • The mathematical tools developed and refined in 
Talasila, Clemente-Gallardo, van der Schaft and Desbrun, Hirani, Leok 
and Marsden j2] could prove useful in developing a more rigorous, coordinate- 
invariant formulation of DTH dynamics and SEM integration. 
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Appendix 
Lemma [B] 

Proof. Using implicit differentiation, we have ||z A || = \\fe l \JH z \\ < M x . 
Therefore ||^(A 2 , ^) — z(Xi, z k ) \\ < M x |A 2 — Ai| . Using the abbreviation f z ( A) 
for f z (X, z k , z(X, z k )) we have 

WhM - h(x x )\\ = l - \\x 2 m zz {z 2 ) - x x m zz {z x )\\ 

< - |A 2 | \\H zz {z 2 ) - H zz &i)\\ + - \\H zz &i)\\ |A 2 - Ai| 

< -XsJh 11*2 - zi\\ + ~M 2 |A 2 - Ai| 
<-(X slH M x + M 2 )\X 2 -X x \. 

Therefore, 

\\z x (X 2 , z k ) - z x (Xi, z k ) || = \ \\f T \X 2 )JH z (z 2 ) - f T \X 1 )JH z (z x )\\ 

< \ ||/f 1 (a 2 )|| \\n z (z 2 ) - n z (zoii + \ \\n z (z x )\\ \\f=\x 2 ) - j=\x x )\ 

< M 2 \\z 2 - z x \\ + l -M x \\fi\X 2 )\\ \\f z {X 2 ) - fc{X x )\\ \\f^{X x )\\ 

< M 2 M X \X 2 - Ai| + M x (X 8lH M x + M 2 ) \X 2 - X x \ 

< (2M X M 2 + Ml) \X 2 - Ai| 
= lz \X 2 - Ai| . 
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Lemma [HI 

Proof. Using Lemma El and the fact that h(X,z k ) = ip(z(X,z k )) we have 
dh(X 2 ,z k ) dh(Xi,z k ) 



dX 



dX 



= \lp z {z 2 )z\{X 2 ,Z h ) -1p z (zi)Zx(Xl,Z k ) 

< \\^J(z 2 )\\ \\z\(X 2 ,Z k ) -Z\(\i,Zk)\\ 

+ \\M^Zk)\\ HvK^)-^^)!! 

< (N llz + M*N 2 ) |A 2 — Ai| 
= lh \X 2 - Ai| . 



Lemma 1101 (Monotonicity) 
Proof. Inequality (jSJ of Lemma |H1 implies 



dg_ 
dX 



+ -A ( ip k + -ip' k X 



< AK |A| 



(20) 



Under the assumptions of claim (i), we have that |^^A| = \\^' k \ |A| < 
\/2AK \ipkW\ipk\ ]96K = l\i>k\- Therefore, ~ \if) k \ > ||^a;A| , and we have 



IV>A 



< 



(21) 



Now, assume dg/d\ = for A ^ 0. Then (EOJ) implies J |A| |^ fc + |^A| < 
4K |A| 3 < (\ip k \ /96K) \X\ = ± \ip k \ \X\ . Therefore, 



+ ^^feA 



<> 



(22) 



Inequalities (j2~Tj) and (|2^|) imply | < |, a contradiction. Therefore dg/dX ^ 
for A and hence g(X,z k ) is monotonic increasing/decreasing on the 
intervals (— A k , 0) and (0,Afc). 

Under the assumptions of claim (ii), inequality pHjl becomes 



dg{s,z k )/dX lip k 



Assume dg(s, z k )/dX = for s 7^ 0. Then becomes 



1 

< — 




- 12 





(23) 



s| < 



(24) 



If s < 0, (J23|) implies 2 — s < — |s or s > 6, a contradiction. If < s < 2, then 
(J24|) implies 2 — s < |s or s > |. If s > 2, then we have s — 2 < |s or s < 6. 
Thus, dg(s, z k )/dX can equal zero only if s = or | < s < 6. Therefore, 
g(s,z k ) is monotonic increasing/decreasing in the interval (— S k ,0), in the 
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interval (0, Sk) if Sk < §, in the interval (0, §) if Sk > | and in the interval 
(6, S k ) if S'jfc > 6. 

Finally, under the assumptions of claim (iii), if dg/dX = for A ^ 0, 
inequality (|2H becomes § |^| A 2 < AK |A| 3 < 4fsT(|^| /A8K)X 2 = ± |^| A 2 
which implies ~ < ^, a contradiction. Therefore dg/dX is nonzero for A 7^ 
and g(A, z^) is monotonic increasing/decreasing in the intervals (—A*., 0) and 
(0,A fc ). ■ 



Theorem [T2l 

Proof. We have from equation © of Lemma 03 that 



^ Wfc 8 24 fa 



- iVfcl ' ' - 48 



fa 



Using the reparametrization A = — (ipk/fak) s we h &ve 



9{s,Zk) 



*Hk \ffa 

fa s VVi 



For — Sfc < s < we have 



s 2 + 



24 \>£ 



< 



48 \fa k 



fa 



fa 



|A| 3 . 



V± _ 1 A 2 (2 _ s ) < 9(s,z k ) < ?4 _ 1 A 2 (6 
16 48 



(25) 



For < s < Sk we have 



-— A 2 (6-s) < 
fa 48 1 '- 



g{s,z k ) . W 



< 



- -A 2 (2 



ip k fa 16 
Now, by inequality (|25Jl . if — 5^ < s < and TCk/ipk < 0, then 



(26) 



g(s,z k ) 

fa 



Similarly, by inequality 



< 



fa 



1 

48 



A 2 (6+ \s\) < 0. 



(27) 



if < s < 2 and H k /4>k < 0, then 



9(s,z k ) 



< 



Hi 



fa 



-^A 2 (2-| S |)<0. 



(28) 



Since g(0,z k )/fa — "Hk/fai inequality (|27jl and (J28j) imply the following: 
If Tlk/fa is strictly less than zero, then no solution exists in the interval 
(—Sk, 2) establishing claim (i). If 7i k /fa equals zero, then s = is the only 
solution in (—Sk,2) establishing claim (ii). 

Next, we use the Intermediate Value Theorem to establish claim (iii). 
Inequality (|2~%j) implies that for < s Q < 2, 



9{s ,z k ) 
fa 



< 



fa 



16 



A 2 (2 - s ) < 0. 



(29) 



20 



Assume S k > 6 and -^A k (6 — Sk) < 7~Ck/4>k < 0. Using inequality (pjj) we 
have 

° < ^ - 48 Afc(6 ~ ^ — 0^" (30) 

Inequality ()29|) and ()30|) and the Intermediate Value Theorem imply there 
must exist a solution s\ G [2, Sk) establishing claim (iii). 

Proceeding in a similar fashion, if < Hk/ipk < ^^-^(6 + S k ) and — Sk < 
s < 0, inequality (j23|) implies 

g(S k ,z k ) < 74_ l 2( 6 + 3)<Qm 

Since g(0,z k )/ip k = H k /ipk > 0, the Intermediate Value Theorem implies 
that there exists a solution G (—S k ,0). The monotonicity of g(s,z k ) by 
Lemma fTTTlf ii) implies the solution is unique in (—Sk, 0) and claim (iv) is 
established. 

Now assume S k < § and < Hk/ipk < je^K 2 ~ Sk). Inequality (l2Tfl) 
implies 

^<^-i^(2-5 fe )<0. 

Since g(0,z k )/ipk = Hk/ipk > 0, there exists a solution which by the 
monotonicity of g(s, z k ) (Lemma^fii)) is unique in (0, S k ) establishing claim 
(v). To establish claim (vi), assume < Ttt/ipk < yfs {^kl^'kf ' ■ H Sk > |, 
by inequality (J2HJ) 

9(%*k) < Hk_ J>_(th) 2 <0 



^ Vfc 125 

which implies there exists a solution s k , which by the monotonicity of g(s, z k ), 
is unique in (0, |) establishing claim (vi)a. Moreover, if Sk > 6, inequality 
implies 

0<^<^-±AU6-S k )< 9(Sk > Zk) 



i>k ~ i>k 48 Ky ~ ^ k 

which implies there must exist another solution s* k G establishing 
claim (vi)b. Finally, the minimum value of Ti, k /ip k — ^A 2 (6 — s) for s > is 

Tik/ipk - | {^k/'ip'k) 2 ■ If ftk/^k > | {^k/^'kf , then using inequality fl2f| we 
have that for all < s < S k , 

TC k 2 / ib k \ 2 TCk 1 , 0/ s q(s,z k ) 
0< — -- -f <— A (6 - s < yv ' . 

Therefore, no solution can exist on (0, Sk), establishing claim (vii). ■ 
Theorem [Till 

Proof. Since |A| < < \tp' k \ /48K, we have from equation © of Lemma El 
that 

< * IA| 4 <ilAI* 



g(\,z k ) (H k _ ]_ x3 



fk Wk 24 



1^1 ' ' 48 
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If A < 0, we have 



If A > 0, we have 



Tj± _ J_^ 3 < g(X,z k ) 



f k 48 



< ^ - 1a 3 . 



16 



(31) 



_ J_, 3 < g{X,z k ) 
V h 16 " tf k 



< 74 _ 1 3 

- ^ 48 



(32) 



To establish claim (i), assume — ^A 3 . < H. k /^' k < 0. Then if A < 0, inequality 
fl3*Tj) implies 

74 



< 



1 A s ^ g(-^k,Zk) 



Since g(0, z k )/ip' k = TL k /ip' k < 0, the Intermediate Value Theorem implies 
there must exist a solution A^T G (— Afc,0). Uniqueness follows from mono- 
tonicity (Lemma HoT m)). Inequality ()32|) implies that for all < A < A k , 



g{\z k ) 



< - 



< 



and thus no solution can exist on (0, A*) establishing claim (i). A parallel 
argument establishes claim (ii). The lower bound of inequality (|31|) and the 
upper bound of inequality ()32|) establishes claim (iii). Finally, to establish 
claim (iv), assume TL k /ip' k < — j^A k . If — A fe < A < 0, inequality (J3"T|) implies 

q(X,z k ) 74 1 ...a 1~L k 1 . o 

< — ^ H A < — r H A? < 
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and no solution can exist on [— A&,0]. Since by assumption, TC k /i/; k < 0, the 
upper bound of inequality implies no solution can exist on [0, Afe]. If, on 
the other hand, Ti, k /ip k > ^A 3 ,, a parallel argument also implies no solution 
can exist and thus claim (iv) is established. ■ 
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